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Abstract 

A first-order phase transition separating the smooth phase from the crumpled one 
is found in a fixed connectivity surface model defined on a disk. The Hamiltonian 
contains the Gaussian term and an intrinsic curvature term. 
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1 Introduction 



There has been a renewal of interest in membrane physics [1] of Helfrich and 
Polyakov-Kleinert [2,3,4]. The phase transition between the smooth phase 
and the crumpled phase is one of the interesting topics in the surface model 
[5,6,7,8,9,10,11,12]. The smooth (crumpled) phase is appeared in the model 
with infinite (zero) bending rigidity b which is included in the Hamiltonian 
S such that S = Si + bS 2 , where Si and S 2 are the Gaussian energy and the 
bending energy respectively. 

We will concentrate on the phantom and tethered surface model in this Let- 
ter. Models on triangulated surfaces can be classified into two groups. One 
of them is the fluid model which is defined on dynamical connectivity sur- 
faces [13,14,15,16,17,18,19,20], and the other is the tethered model on fixed 



Email address: koibuchi@mech.ibaraki-ct.ac.jp (Hiroshi Koibuchi). 



Preprint submitted to Elsevier Science 



2 February 2008 



connectivity surfaces [20,21,22,23,24,25,26,27,28,29,30]. There is another clas- 
sification of surface models; a model is called real or phantom according to 
whether the surface is self- avoiding or not. The phantom model is still inter- 
esting because the self-avoiding interaction is considered to be irrelevant in 
the limit N^oo [31]. 

We consider that the phase structure of the tethered surface model remains 
to be studied further. It has been recognized that the phantom tethered model 
undergoes a second-order phase transition both on a sphere [20,21,22,23,24,25,26,27,28] 
and on a disk [29,30]. On the other hand, it has also been reported that there 
is a first-order transition in the model with Hamiltonian slightly different from 
the ordinary one of Helfrich and Polyakov-Kleinert on a sphere [25] . First-order 
transitions can also be seen in a model of Nambu-Goto Hamiltonian with a 
deficit angle term [26] and in a model with Hamiltonian containing the Gaus- 
sian term and an intrinsic curvature term [20]. Therefore, phase transitions of 
the tethered surface model should be understood more clearly. 

In this article, we would like to investigate the dependence of the phase tran- 
sition on the existence of boundary in the surface. For this purpose, we define 
the model on a triangulated disk, which is topologically identical with a sphere 
with a hole, The disk is obtained by dividing the hexagon. The Hamiltonian is 
given by a linear combination of the Gaussian term and an intrinsic curvature 
term; intrinsic curvatures in membranes have no direct analog in one dimen- 
sional objects like linear polymer as noted in [1]. It was reported that the 
model with intrinsic curvature undergoes a first-order transition on a sphere 
[20] as mentioned above. Hence, it is interesting to see whether the boundary 
influences the first-order phase transition of the model. 



2 The model 



The model is defined on a triangulated disk which is obtained by dividing the 
hexagon. Figure 1 is a lattice of size N — 61, which is obtained by dividing 
every edge of the original hexagon into 4-pieces. If the edges are divided by 
L-pieces we have a lattice of (N, N E , N T ) = (3L 2 +3L+1, 9L 2 +3L, 6L 2 ), where 
N, Ne, Nt are the total number of vertices, the total number of edges, and 
the total number of triangles. 

The Gaussian tethering potential Si and the intrinsic curvature S3 are defined 
by 

s^^Xi-Xj) 2 , s 3 = £(^-<5 ) 2 , (i) 

(ij) * 
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Fig. 1. A triangulated disk of size (N,N E ,N T ) = (61,156,96), where N,N E ,N T 
are the total number of vertices, the total number of edges, the total number of 
triangles. 

where J2(ij) is the sum over bond (ij) connecting the vertices Xi and Xj, and 
5i in S3 is the sum of angles of the triangles meeting at the vertex %. The value 
of the symbol 5o in S3 is given by 

2ir, (internal vertices); 

= { (2) 
7r(or 2n/3), (boundary vertices). 



27r/3 is assigned to the six vertices, and 7r to the remaining vertices on the 
boundary. The definition of 5 in (2) allows us to call S3 as a deficit angle 
term. 

The partition function is defined by 



N 

Z(a)= / n^exp[-SpT)], (3) 
S(X) = S 1 + aS 3 , 

where N is the total number of vertices as described above. The expression 
S(X) shows that S explicitly depends on the variable X. The coefficient a is 
a modulus of the intrinsic curvature. The surfaces are allowed to self- intersect. 
The center of surface is fixed in Z(a) to remove the translational zero mode. 

It should also be noted that the Gaussian tethering potential Si at the internal 
vertices differs from the one at the boundary vertices because the co-ordination 
number of them is different from each other. The intrinsic curvature S3 at the 
internal vertices can also be different from those at the boundary one on 
fluctuating surfaces due to the same reason for Si. These differences of Si and 
S3 between the internal and the boundary vertices are the one we can not find 
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in compact surfaces such as a sphere. 

We comment on a relation of S3 in (1) to an intrinsic curvature energy 
S3 = - Ei log(5i/2vr) in [20]. This S 3 in [20] is intimately related with the 
co-ordination dependent term S 3 = — Ei l°g(?i/6), which comes from the in- 
tegration measure Fli dXiq® [32] in the partition function for the model on a 
sphere. S3 = — Ei log(<5i/27r) becomes minimum if Si = 2ir for all i, and hence 
it becomes smaller on a smooth sphere than on a crumpled one. However, the 
flat configuration of the hexagonal lattice shown in Fig. 1 does not minimize 
this S 3 , because there are vertices such that Si ^ 2ir. This is the reason why 
we use S3 in (1), which becomes zero on such flat configuration as shown in 
Fig. 1. 



3 Monte Carlo technique 

The canonical Monte Carlo technique is used to update the variables X so 
that X' = X + SX , where the small change SX is made at random in a small 
sphere in R 3 . The radius Sr of the small sphere is chosen at the beginning 
of the simulations to maintain the rate of acceptance rx for the X-update as 
0.4 <r x < 0.6. 

We use a random number called Mersenne Twister [33] in the MC simulations. 
Two sequences of random number are used; one for 3-dimensional move of 
vertices X and the other for the Metropolis accept/reject in the update of X. 

We use surfaces of size iV = 2107, iV = 3997, iV = 6487, and iV = 9577, which 
correspond to the partitions L = 26, L = 36, and L = 56 respectively. The size 
of the surfaces is relatively larger than that of spherical surfaces used in [20]. 

The total number of MCS after the thermalization MCS is about 1.8 x 10 8 at 
the transition point of surfaces of size iV = 9577, iV = 6487, and iV = 3997, and 
1.6xl0 8 for iV = 2107. Relatively smaller number of MCS (1.0x10 s ~ 1.6xl0 8 ) 
is done at a distant from the transition point. 



4 Results 

Scale invariance of the partition function predicts that Si/N — 1.5, which is not 
influenced by the boundary condition imposed on the model. We also expect 
that this relation should not be influenced by whether the phase transition is 
of first order or not. Figures 2 (a) and 2 (b) are Si/N against a obtained on 
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surfaces of size N = 2107, N = 3997 and N = 6487, N = 9577. We find from 
these figures that the expected relation Si/N—1.5 is satisfied. 
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Fig. 2. Si/N vs. a obtained on surfaces of (a) iV = 2107, iV = 3997, and (b) iV = 6487, 
iV = 9577. 



We expect that the surfaces become flat at a — > oo and crumpled at a — > 0. 
Hence the size of the surfaces can be reflected in the mean square size X 2 
defined by 



X' 



(4) 



In order to see the size of surfaces, we plot X 2 against a in Fig. 3 (a). We see 
that X 2 disco ntinuously changes at finite a. This indicates a discontinuous 
transition between the smooth phase and the crumpled one. 
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Fig. 3. (a) X 2 vs. a obtained on surfaces of size from N = 2107 to N = 9577, and 
(b) log-log plots of X 2 vs. N obtained just below and above the transition point a c 
on each surface. The straight lines in (b) are drawn by fitting the data to Eq. (5). 
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Figure 3 (b) shows log-log plots of X 2 against N, obtained just below and 
above the transition point a c where X 2 discontinuously changes as shown in 
Fig. 3 (a). The straight line denoted by smooth (crumpled) is obtained by 
fitting X 2 in the smooth (crumpled) phase just above (below) a c by 

X 2 ~ N 2 / H , (5) 

where H is the Hausdorff dimension. Thus we have 

f 1.86 ±0.09, (smooth); 
H={ (6) 

[22.1 ±13.2, (crumpled). 



The result H = 1.86(9) in (6) implies that the surface is almost smooth in 
the smooth phase. On the contrary, H = 22.1(132) implies that the surface is 
highly crumpled in the crumpled phase. 
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Fig. 4. Auto-correlation coefficient A(X 2 ) in the smooth and the crumpled phases 
on surfaces of size (a) iV = 3997, and (a) iV = 9577. 

The phase transition can be reflected in the convergence speed of the MC sim- 
ulations. In order to see this, we plot in Figs. 4(a) and 4(b) the autocorrelation 
coefficient A(X 2 ) of X 2 defined by 



MX-) = S*W*»>, (7) 
Ei* 2 (n)] 2 

Tj+i = Ti ± n x 500, n = 1, 2, • • • . 

The horizontal axes in the figure represent 500xn {n = 1, 2, • • -)-MCSs, which is 
a sampling-sweeps between the samples X 2 {ti) and X 2 (r i+ i) . The coefficients 
A(X 2 ) in Figs. 4(a) and 4(b) are obtained on surfaces of size iV = 3997 and 
N = 9577 respectively. We clearly see from the figures that the convergence 



6 



speed in the smooth phase is about 10 times larger than that in the crumpled 
phase. The reason why the convergence speed in the smooth phase is so larger 
than in the crumpled phase is that the phase space volume (C R 3 ), where the 
vertices X,- L take their values, is relatively larger in the smooth phase than that 
in the crumpled phase. This difference in the phase space volume can also be 
reflected in X 2 shown in Figs. 3(a) and 4(b). 

However, we find no phenomena of the critical slowing down in A(X 2 ) at the 
transition point. The reason of this is because of irreversibility of the transi- 
tion. The change from the crumpled phase to the smooth one appears to be 
irreversible. If the vertices once localize to a small region in R 3 at the tran- 
sition point, they hardly expand to be a smooth surface. This irreversibility 
was also seen in the same model on a sphere [20]. 
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Fig. 5. (a) S3/N vs. a, and (b) S2/NB vs. a, where Nb{=Ne~ 6L) is the total 
number of internal bonds. 

Figure 5(a) shows the intrinsic curvature energy S 3 /N changing discontinu- 
ously against a. This indicates the existence of a first-order phase transition. 
The variation of S^/iVe against a is shown in Fig. 5(b), where S2 is the bend- 
ing energy defined by 62 = I]j(l — cos 9i) and Nb is the total number of internal 
bonds. It should be noted that 6L-bonds (the total number of bonds on the 
boundary) are not included in N B (=N E — 6L). The quantity 1 — cos#j is de- 
fined only on the internal bonds. The discontinuous change of S2/NB shown 
in Fig. 5 (b) indicates that the phase transition separates the smooth phase 
from the crumpled phase. 

Snapshots of surfaces are shown in Figs. 6(a) and 6(b), which are obtained at 
a — 13 and a— 13.5. The size is iV = 9577 in Figs. 6(a) and 6(b). Sections of the 
surfaces in Figs. 6(a) and 6(b) are shown in Figs. 6(c) and 6(d) respectively. 
It is easy to see that the surface in Fig. 6(a) is crumpled and that the surface 
in Fig. 6(b) is smooth. Although the surface in Fig. 6(b) looks like a cup, it is 
almost smooth as can be seen from the section in Fig. 6(d). 
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Fig. 6. Snapshots of surfaces of size N = 9577 at (a) a = 13, and at (b) a = 13.5, 
and (c) a section of the surface in (a), (d) a section of the surface in (b). Figures 
are shown in the same scale. 

5 Summary and conclusion 

We have shown that a tethered surface model with an intrinsic curvature 
undergoes a first-order phase transition between the smooth phase and the 
crumpled phase on a disk. Our results indicate that there is no influence of 
the boundary condition on the phase transition of the model with intrinsic 
curvature. 

The Gaussian tethering potential and the intrinsic curvature energy are de- 
fined on the triangulated disk which is obtained by dividing the edges of the 
hexagon into L-pieces. The triangulated disk is characterized by (N, N E , N T ) = 
(3L 2 +3L+1, 9L 2 +3L, 6L 2 ), where N E , N T are the total number of edges and 
the total number of triangles. We used surfaces of size up to N = 9577, which 
is about four times larger than those used in [20] for the similar model on a 
sphere. Hence, the observed first-order transition in this Letter is considered 
as the one in the thermodynamic limit of the model. The canonical Monte 
Carlo technique was used to update the position variable X of vertices. 

Combining the result in this Letter and those in [20], we can conclude that 
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1) the tethered surface model with intrinsic curvature undergoes a first-order 
phase transition 

2) the phase transition can be observed both on a sphere and on a disk. 

It is also interesting to study the phase structure of the tethered model on a 
torus, and that of the fluid model. 

This work is supported in part by a Grant-in- Aid for Scientific Research, No. 
15560160. 
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